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Abstract 

In this note, we derive the spectrum of the infinite quantum XXZ spin chain with domain 
wall boundary conditions. The eigenstates are constructed as limits of Bethe states for the 
finite XXZ spin chain with li^ (sk) invariant boundary conditions. 
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1 Introduction 

In this paper we consider the XXZ spin chain with domain wall (DW) boundary conditions. 
In physical terms it describes a one-dimensional lattice where each point carries a spin in- 
teracting with its neighbors. Its time evolution is given by the Hamiltonian 

H = -\ E KV+i + ^+i + ^ir 4 )' (i-i) 

feZ+l/2 

where k denotes the position on the lattice, oi is the z-th Pauli matrix acting on the spin at 
position k, and A is the parameter characterizing the anisotropy of the spin interaction. In 
our case, A > 1. 

This is one of the most studied systems in statistical physics. It has a number of remark- 
able properties. In particular it is integrable, i.e. the Hamiltonian ( |1.1| can be diagonalized 
using the Bethe ansatz [JJ. 

In this paper we study the spectrum of xxz Hamiltonian in the space of states with do- 
main wall boundary conditions. This space is the natural zVcompletion of 

span { f\ a~ Y\ o-+Cl DW | x t < 0, y f > } , (1.2) 

i i 

where Q D w = |- • • Ttil ■ ■ ')■ 

This paper is the continuation of [2j, where it was shown that the evolution of a system 
of random integer partitions can be mapped to the half-filled sector of the ferromagnetic 
XXZ spin chain with domain wall (or kink) boundary conditions. In [2], the quantization of 
the system of integer partitions or Young diagrams was considered. There is a well-known 
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Figure 1: Two-dimensional partition corresponding to the state flow ( [10, 6, 3, 3, 3, 2, 1, 1] ) = 

V W V 7 - I5tp* 9 ip_9 Xpt tp_ 5 Q DW . 

2 2 2 2 2 2 

map from partitions to the NS sector of a one-dimensional free fermion. A fermionic creation 
operator tpl corresponds to each black dot of the Maya diagram, i.e. the projection on the hori- 
zontal line (see Fig.[l). The empty partition Qdw is obtained by filling all negative positions. 
Any other partition can be obtained by acting with creator-annihilator pairs on Qdw- hi a 
quantization procedure familiar from the theory of quantum dimers (see EHU), the Hilbert 
space of the quantum system is spanned by vectors that are in one-to-one correspondence 
to the classical configurations (the integer partitions). In terms of spin chains, this Hilbert 
space is precisely J^bw- The ground state Boo of the quantum Hamiltonian is required to re- 
produce the steady state distribution of the classical system. A natural Hamiltonian fulfilling 
this requirement is given by 

\ 

H = - $t+i#fc + fk^k+i ~qn k (l- njfc+i) - - tt fc+ i (1 - n k ) , (1.3) 

icez+i/2 9 

where n k = ip? tp k is the fermion number operator. This Hamiltonian can be recast into the 
form of Eq. ( |1.1[ ) by using the Jordan-Wigner transformation. The quantization of the growth 
of integer partitions thus corresponds to the XXZ spin chain with domain wall boundary 
conditions. 

The quantum xxz model is also related to a classical random process system, to the asym- 
metric exclusion process (asep). This is a Markov process on a one dimensional lattice with 
particles on its sites. The particles diffuse with different probability rates for jumping left or 
right. None of the sites in such a system can be occupied by more than one particle. The tran- 
sition rate matrix W for such a process can be obtained by conjugating the xxz Hamiltonian 
above with a diagonal matrix. 

The gap in the excitation spectrum has been determined in [5j. We construct the whole 
spectrum. Because the transition rate matrix is symmetrizable, the ASEP is invertible and 
therefore it converges to the equilibrium state which is the ground state of the xxz Hamilto- 
nian over the DW space of states. 

The problem of computing matrix elements (form factors) 

(1-4) 

for all eigenvectors xp a , xpp of H remains. 
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Our main results are the following: we describe the eigenvectors and eigenvalues of the 
XXZ Hamiltonian on the DW Hilbert space and we find explicit expressions for the form 
factors of the ground state extending results of [2J. We also study the q — > 1~ limit of these 
(scaling) matrix elements. 

The plan of this paper is as follows. In Section |2j the recall some basic facts about xxz 
Hamiltonian, and the asymmetric exclusion process (ASEP) is introduced. The main part of 
the paper is Section]3] where we describe the spectrum of the xxz model with DW boundary 
conditions. In Sec. [4} the time evolution for the ASEP is treated. In Sec. [5} some calculations 
are performed on the ground state, such as the limit shape and some matrix elements. 



2 The xxz Hamiltonian 

Here, we recall the definitions of the DW Hilbert space and the ASEP Markov process. 



2.1 The DW Hilbert space 

Denote by Q the ferromagnetic ground state with all spins up, and by Cl(x\, . . . ,x m ) the state 
where the spins at positions X\ < Xj < . . . are down. 

The domain wall (DW) ground state flow is the state where all the spins with negative 
coordinates are up and all the spins with positive coordinates are down. The states 

n m 

n^n< Q ow (2.i) 

1=1 z=l 

form a basis in the DW Hilbert space. When n = m, these vectors can be parametrized by 
partitions as Dpiy where 3 = \5\, . . . , S r ] S is the integer partition which has u and v 
as modified Frobenius coordinates: 

Vi = 5 i -i+ 1 - m = - U - i + ^ , (2.2) 

where i ranges from 1 to the number m of the squares on the diagonal of the Young diagram 
of 3. The partition 3* is the reflection of 3 with respect to the diagonal. 



2.2 The ASEP evolution of the DW states 
2.2.1 Continuous Markov process 

Recall that a continuous-time Markov process is defined by the following evolution equa- 
tion: 

^l=Y f W(C,C')Pc(t), (2.3) 

where C is a configuration of the system, Pc(t) is the probability that the system is in config- 
uration C at time t, and W(C, C) is the transition rate from C' to C. Because of yj c Pc(t) = 1/ 
and that Pc(t) > 0/ one finds that W(C, C) must satisfy the following conditions: 

W(C,C')>0 VC'/C (positivity), (2.4) 

W{C,C) = ~ W{C',C) (total probability). (2.5) 
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Here we will focus on Markov processes satisfying the detailed balance condition (so-called 
invertible Markov processes). The transition rates in such processes satisfy the symmetry 
relation 

W(C', C)P£ = W(C, C')P° C , (2.6) 

where P® > and we choose them to be normalized as P® = 1. The symmetry of the 
transition matrix implies that {P°} is an eigenvector of W with eigenvalue 0: 

Y^W(C,C')P°c =0, (2.7) 

c 

One of the fundamental features of invertible Markov processes is that such a process start- 
ing from any distribution P(0) converges to P° as t — » oo. 

2.2.2 Asymmetric exclusion process 

An asymmetric exclusion process (ASEP) is a Markov process on the infinite spin chain. In 
this paper we will focus on the study of such processes in the space J^bw- 

The transition rate matrix for the ASEP is determined by its off-diagonal matrix elements 
which can be described loosely as 

w({... Tl... },{... IT ...}) = or, w({... IT... },{... Tl...}) = o> + . (2.8) 

This linear operator can be written as 

W = £ (oT a+a k+1 + w+ + H^+HL (afaf +1 - l) V (2.9) 



Let us parametrize the weights w ± by w =t = A ^ ±x . Then it is clear that the transition rate 
matrix for the ASEP has the symmetry 

W(A, q- 1 ) = CW{A, q)C , (2.10) 

where 

c =®(! J)- p- 11 ) 

Because of this symmetry we will assume q < 1 for the rest of the paper. 
Define the linear operator 

£))•(§(*!))■ 



!<0 v ' 7 / \i>0 



It acts on the space J^bw as 



m n m n 



UU^IKPdw = q^-^Yl^lla+noy,, (2.13) 

i=l ;'=1 i=l /=1 

and it is a bounded operator acting on this Hilbert space. It is easy to check that 

UafU' 1 = q M vf . (2.14) 
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These identities have two important implications: 

W = -AUHU' 1 , W = U 2 WU- 2 , (2.15) 
where W f is the transposed matrix W and 

« = " E (Vfc+i + + f " 1)) • (2-16) 



Here A = + ^ x ). Therefore, the matrix W is conjugate to the XXZ Hamiltonian and it 
satisfies the detailed balance condition. This implies that the vector 

= ^ E E 9 OU<w-*> IK n<«D W (2-17) 

)!>0 xi<ar2<...x«<0<yi<-"<i/,! z=l ;'=1 

is a normalized eigenvector of W with eigenvalue and that the ASEP evolution with the 
transition rate W converges to Yo for any initial condition as t — > oo. 

3 The Bethe ansatz 

In this section, we construct the full spectrum of the Hamiltonian on the Uq(sl2) -symmetric 
XXZ chain. 

3.1 The Bethe ansatz for the XXZ chain over the ferromagnetic vacuum 

The Bethe ansatz [lj for the xxz spin chain describes explicitly the eigenvectors and eigen- 



values of the Hamiltonian in Eq. (1.1 1 over the ferromagnetic vacuum Q. More specifically, 
consider the vector 

B(z)= £ A(x|z)Q(x), (3.1) 

Xi< — <X m 

where 

A(x|z)= £ e n A(z nil) ,...,z nim) )z x ^ iy ...,z x ^ m) = £ e n A(z n )z x n . (3.2) 

Here |z ; | = 1, the sum runs over the permutations of {z\,...,z m }, and e n is the sign of n. 
The explicit expression for A(z) is 

A( Zl ,...,z m ) = YlA( Zi , Zj ) = Yl(l-2Azi+ZiZj) . (3.3) 



Bethe proved that this vector is an eigenvector of ( 1.1 1 with the eigenvalue 



£(z) = £(2A-z a -z- 1 ) . (3.4) 



a=\ 



Such eigenvectors are called m-magnon states. 

In sectors with m > 2, the magnons can form bound states. These bound states decay 
exponentially with the distance between the magnons. In the eigenvector corresponding to 
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a bound state of n magnons, the numbers {z\, . . . , z„) G (C*) n are such that \z\ . . . z n \ = 1, 
\zj\ > \zi + i\, and Z\ = h n ~ l (z n ), where h(z n ) = 2A — The corresponding eigenvector is 

x r x i 

x^-\ — \-xn / Zi \ n 

B(zx,...,z n ) = E A{zx,...,z n ) (zi...z„) «~ Yl ( ) Cl{x lr ...,x n ) . 

x 1 <x 2 <---<x„ \<i<j<n V z !/ 

It is convenient to use the parametrization of Z; by rapidities, see e.g. |6l: 



(3.5) 



_ sinh(f/3 + (f-; + l) 7 ) _ sinh(fg + f) 

Z '~ sinh^ + (§-;)?/) ' Zx ~ n_ srnh(z/3-f)' 1 J 

where A = cosh(//). In this parametrization, the energy of a bound state of n magnons is 

_ 2sinh(?y) sinh(n?y) 
MW ~cosh(m;)-cos(2fi)- ^ 

The eigenvectors with an arbitrary number of magnons and their bound states form the 
complete continuous spectrum of the XXZ spin chain over the ferromagnetic vacuum. 

3.2 The spectrum of the XXZ Hamiltonian over the DW vacuum 

The Bethe states above do not belong to the Hilbert space ^bw- It is however possible to 
map them to a complete set of solutions of the XXZ Hamiltonian for the infinite chain with 
domain wall boundary conditions. 

Our main result is the following. The states Boo(z) defined by 

Boo(z) =q^> ?' 1 E r™ & +7 q dw(<$) (3.8) 

form a complete set of eigenstates of the continuum spectrum for the XXZ Hamiltonian for 
the infinite chain with domain wall boundary conditions with eigenvalues 

m 

f(z) = X:(2A-z fl -2- 1 ) . (3.9) 

a=i 

Here, the z are as above, 5? is the set of integer partitions, \8\ = 8\ + • • • + 8 r , and the A(x|z) 
are defined in Eq. ( |3.2| . 

The immediate corollary of this is that the vectors 

o,oo = ub„(z) = e ^ E A( ^" ?1 ;:-/;r ?m|z) Q ^w (3.io) 

form a complete set of eigenstates for the rate matrix W = AUTiU^ 1 of the ASEP process. 

The strategy of the proof is to pass to the limit N — > oo from the spectrum of the XXZ 
Hamiltonian for finite spin chain with XlAslz) invariant boundary condition. The latter was 
computed in [7]. Let (z) be a Bethe vector for the spin chain of length N (see Appendix 
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B). The Hamiltonian is invariant with respect to the diagonal action of U^slj), therefore the 
vectors 

( p ) 

B[ N) (z) = ±^K- p/2 FPBW (z) , (3.11) 

[P\q 2 - 

are also eigenvectors of such spin chain. Here 



1 Cj 



and K and F are the diagonal actions of the generators of U^sh) (see appendix A) on C 20N : 

K = q°* ® • • • (8) ^ , (3.13) 

F = £]^ 3 (8) • • • ® ^ ® a~ ® 1 (8) • • • ® 1 = £]L y cr- . (3.14) 

y " r; ' y 

Lemma. We /iai?e the following identity: 

^® T K- p/2 FPQ(x 1/ . . . ,* m ) = ^ (?P (2 P +2m-l)+Ef =1 (i/ ! -2/( !/i |x)) x 

xO({^}^{y ; }P), (3.15) 

w/jere ? (y z - 1 x) = #{x a \x a < J/;}. 

Proof. The commutation relations between a~ and Ly are given by 



y L ?/2% ifyi<y2- 



^= t w i y 1 ^ 1 (3^) 



From here we find: 



F v — L yi o- yi L y2 o~y 2 . . . Ly p o-y v — q 2 '^ L Vl Ly 2 . . . Ly p cr yi . . . a yp 

= iP}q 2 ! E L Vi L l/ 2 • • • L y P % • • • V 

yi>y 2 ,->y v 

where n is a permutation in the symmetric group S p , and /(zr) is the length of the permuta- 
tion. With this, 

x m ) = [p]^! X] L yi ...L yp O({x fl }^ { Vl }\), (3.18) 

yi>->y p 
yi^« 

where {x a }™ denotes the positions of the m particles of the original state and {yi}{ denotes 
the set of positions of the spins flipped by the action of the cr - . The lemma follows after 
evaluating the action of L yr □ 
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It is convenient to introduce new coordinates. Assume that y,- < J/j+i, Xj < Xj+i, yu-\ < 

x\ < yi v and i/; lH v\^-\ < x a < Vi v --+i a - Write the union of sequences {xi}f and {i/;}^ as a 

sequence p with x fl = U; a/ fl = \,...,m and U/ a+)t = y h ,.. + j a+k ,k = l,...,l a+1 - 1. In 

other words: 

u = {y\, . . . ,xi, . . . ,x 2l . . . ,x m , . . . ,y v } = {ui . . . ,ii h , . . . ,u h , . . . ,ui m , . . . ,u p+m } . (3.19) 



It is easy to check that 



£%ix) = -r 2 +1 ) + E^ 



(3.20) 



i=l 



8=1 



In terms of {«} the formula in Eq. ( 3.17) is 



[p]^ ! 



X 



X 



^ « ! -Er= 1 (%+2/„) Q ( Ml/ . . . , Mp+m ) . ( 3.21) 



N 1 \ \ Nil 



Now we want to pass to the limit N — > oo, p = y — m with finite m. In terms of variables 



Si = Ui-\ + i, 



(3.22) 



the formula for (z) reads 

B (N) ( , _q^_ K - P/2pP 



">Sl>->Sp+m>0 



E 



A( N )(x|z)Q(x) 

\W(Si 1 -h + l,...,6t m -l m + l\z) 



l<h<--<l,n<P+m 

xCl(St-l r ...,S v+m -p-m + l). (3.23) 

Now, set p = y — m. In this case, the summation in the formula above is taken over integer 
partitions 5 G J^n corresponding to Young diagrams which can be inscribed into a N/2 x 
N/2 square. We can rewrite this formula as 



B (N) fzl 



E 



q ("T) + \S\ 



E 



\.W(5 h -h + l S lm -l m + l\z) 



l</l<-<«m<4 



(3.24) 

where |£| = + • • • + S r and Qdw(^) = 0(£i — \,.. ■ ,$n/2 ~ y + \)- Taking the limit 
N — ► oo in this expression is premature. In order to get a convergent series we have to 
rearrange this sum. 

First consider the one-particle Bethe state: 



B (N) (z) _ r ^fA^-l+M 



(3.25) 
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This expression can be written as: 



B 



= E t m 



L, A(") - / + I | z ) n/2 A< N > (-/ + i |z) 



E 

/=1 



+ E 

l=r+l 



Cl DW (5), (3.26) 



where we can use the explicit form of A^ N ) {x\z) (see App.jBj) to evaluate the second sum: 

(3.27) 



n/2 A (N )( _ ?+ i| z) _ lAW(-r + ||z)- A(«) (-r - § |z) 

E - 

;=r+i */ 



Now we can then take the limit N — > oo and and we obtain the formula (3.8 1 for m = 1: 



U\(Si-l\z) - A(-r-l|z) 



E 



+ 



<f + 



^Dw(<5) / 



(3.28) 



Here 5? is the set of integer partitions. Using the explicit expression for A(x|z) = z x we can 
further rewrite Boo as 



Beo(z) = £ q 



£A(^-/|z) (z 9 ) 



Z=l 



1 — qz 



(3.29) 



The identity 



E 



A(^-/|z) 



q s l+ i 



E 



A(^-/|z) 



z=i ■» ' i=i 1 



S,+l 



+ E r^ = E 

Z=r+1 (Z<?J ;=1 



A(^-/|z) (z<?)- 



(3.30) 



which holds for q > 1 demonstrates that the formula (3.291 can be obtained from ( 3.24} by 
taking the limit in the coefficient functions for q > 1 and then by continuing to q < 1 after 
rearranging the sum. 

The same reasoning can be repeated in the general case m > 1 and leads to the following 
expression for Boo(z): 



b-w = q w E ^ E 

sey i<k<...l m <oo 



^Dw(^) / 



(3.31) 



where for q < 1 each divergent geometric series is to be understood in the sense of Eq. {330) . 
3.3 Scalar products 

The scalar products between two eigenstates can be evaluated by considering the commuta- 
tion relations of the Uqfah) algebra (see Appendix [A|. The scalar product is given by 



(B.CABcW) = (B(z'),B(z))n T — » 



z.e. a universal factor times the scalar product of the two primary Bethe states. 



(3.32) 
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To prove this identity consider a chain of length N. Let us evaluate the scalar product 

(B;, N) (z'),l4 N) (z)) = f (r^VBW^r^W). (3.33) 

P [P\q 2 -[P lq 2 - 

First notice that 



and 



The elements £ and F satisfy a well-known identity first derived by V. Kac: 

e p' f p = f p-p> K -p' q y' fl 1 " K2 r 2p+2k 



(3.34) 
(3.35) 



fc=l 



1-q 



+ terms proportional to powers of E . 

(3.36) 

Here p > p'.lfp < p', every term is proportional to some power of E (from the left). 

By construction, the Bethe states are highest weight states and EB(z) = 0. Therefore all 
terms proportional to E with vanish in the scalar product and we obtain 



B 



W ( z '), k-p n f\ 7 B (N) (z) ) 



In the special case N = 2 (p + m) the expression simplifies and reads: 



(B p (p+m) (z'),B^ +m) (z)) = V( B2(P+m) ( z ')^ 2( P +m) (z)) 



2p' 

P. 



(3.37) 



(3.38) 



where [ 2 ^] 2 is the ^-binomial coefficient. We are now in the position of taking the p — > oo 



limit and find: 



(Boo(z')> Boo 00) = lim(B(z / ) / B(z)) 



2p 

P. 



(B(z'),B(z))n T 



(3.39) 



4 Time evolution for the ASEP 



Having obtained the complete set of eigenstates for the XXZ chain with domain wall bound- 
ary conditions puts us in the position to clearly state (if not solve explicitly) the problem of 
the time evolution of the asymmetric exclusion process. Consider a system that at initial 
time is in the configuration Q(°) = Qdw/ where all the spins at negative positions are up 
and the others are down (empty partition state), which evolves with transition rate W. After 
a time t, the system will be in the state n^ W/ given by 

Q^ w = e fW Q DW . (4.1) 

In Sec. [4] we found that W is related to the XXZ Hamiltonian by the similarity transformation 
W = —AUT-CU^ 1 . This means that the time evolution operator can be written as 

e m = Ue~ tAH U-\ (4.2) 
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The exponential of TC can be decomposed in terms of the eigenvectors given in Sec. 3.2 using 
the spectral decomposition theorem: 



n t DW = Yl(i- q 2k )u 



/c=l 



m tW ||B(z)|| 2 

f / dz (WW rt ^ )0co(z) 

tJ \\B(z)f 



n (i-? 

^=1 



, (4.3) 



where £(z) is the eigenvalue ( |3.9| corresponding to Boo(z), and Geo = Boo is the ASEP ground 
state. 



5 The limit shape and matrix elements for the ground state 

The problem of calculating the magnetization profile of the ground state Boo(0) is closely 
related to the classical problem of finding the limit shape of the ensemble of random parti- 
tions, where the limit of vanishing lattice spacing and q — > l - is taken. The limit shape is 
simply the integral of the kink shaped magnetization profile of the ground state. Here, we 
recall results of O, making use of the techniques developed in the present article. 

The magnetization profile of a state of the spin chain labeled by a vector v is given by 

Consider a chain of length N. The ground state fulfills Bj^(O) <x (S~) p Q( N ), thus the corre- 
sponding magnetization profile reads 

(K-^FPbW(z>),o3 x K-^F?bW(z)) 

m B< N) (0)W (K-P/2 F p B (N)( z /) /K - P /2 F p B (N)( z )) • ^ 

This can be evaluated explicitly using the commutation relation 

WlF P ] =2[?V L * E hi--- Vi<V" (5-3) 

y\<--<y v -\ 

and using the fact that = Q^. 

The final expression for N = 2p is given by 



m (0 Jx) = l-2V a 2k{x+ ^ ^ P ' q h 



(5.4) 



It is easy to obtain the limit shape, when we consider first the limit p — > oo for fixed q, and 
then q — > 1 . Sending p — > oo we obtain 

m Ba> (x)=l-2£(-l) k q 2 ®q 2kx . (5.5) 

k=0 
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Assume that the variable u = —2x log q remains finite when q — ► 1 . Then the magnetization 
profile converges to 

CO -1 u 

mB >) = l-2£(-l)V^< = i-^, (5.6) 

and the corresponding limit shape is the antiderivative of — niB^{u) : 

p(u)=21og(2cosh(D) . (5.7) 

Indeed, ]i(u)' = —m^iu), and \i(u) —> \u\ as \u\ — * oo. 

The detailed probabilistic analysis of the limit p — > oo with being fixed is done in l8l . 



Matrix elements for the (projector on) the ground state 

Here we calculate matrix elements for the projector on the ground state, i.e. the probability 
on the ground state to have the Zi-th particle move by 5\, the /2-th particle move by 82, etc. 
starting from flow- 
Consider the half full sector for a chain of length 2p. The ground state can be written as 



B 



(2p) 



E q^mW( yi y p ). 

-p+\<y P <-<y\<p-\ 



(5.8) 



The probability for the Z-th particle from the right to move by 5 is given by 



P{1,5) 



(B^\s( yi -x)B^) 



where x = 5 — l + ^is the position in the chain of the Z-th particle. Explicitly, 



(5.9) 



P(l,S) 



£ 6( yi -x)q 2 ^ 
-p+\<y P <---<y\<p-\ 

E 7^' 

-p+\<y P <--<yi<p-\ 



-■p+l<y p <--<yt + i<x-l 



qlZli+iVi 



E 



x+l<yi-i<-<yi<j>-j 



E 

-p+\<y P <-<y\<p-\ 



(5.10) 



Using the identity 



the probability is 



E </ 

a+l<y m <---<yi<b 



2UUy, 



b — a 
m 



P{U) = q 



218 1 



rp+5-h r-p—S+l— li 

Jq 2 ^ l-l if 



\ 2V \ 



(5.11) 



(5.12) 



12 



and in the infinite case we find 



P(l,S) = q 



216 



(<?W);-i(<? 2 ;<? 2 )<5 



(5.13) 



Note that the magnetization profile in Eq. \5A) can be obtained starting with these probabil- 
ities and summing for / G IN, at fixed position x. 

In a similar fashion one can calculate the joint probability of having the Zi-th particle 
move by Si, the /2-th particle move by S2 and so on. It is given by 



It can be expressed as 



-p+\<y v <--<yi<p-\ L«=i 



-p+\<y P <-<yi<p-\ 



(5.14) 



P({/i}r =1 4^}ti) = ^ l/l+2E&1<5i+l(/i+1 - Z,) 



rp+Si-h rp-<5„+/„-l-| 



[2 P 1 



fl-1 

><n 

1=1 



^ - ^+1 + U+i - k - 1 



, (5.15) 



and in the p — > 00 limit, 



P(l /• 1" f 1 - fl 2 ' J i / i+ 2 ^=i 1 ' 5 m Cm-'/) ; ^ )°° x 

1 U l i / =1 ' l °t fi=l) — H ( n 2. n 2\ („2-„2\ 



n-l 

1=1 



?>i - + U+i - k - 1 



(5.16) 



Note that this is the same as the probability of having an integer partition with the number 
Si at position ?i, the number S2 at position I2, and so on (here it is obvious that Z; < Z i+ i and 

h h In 

P([...Si---S 2 ---S n ...}). (5.17) 
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A The algebra l^(sZ 2 ) 

The algebra UJsh) is generated by E, F, K, and K^ 1 under the relations 

Kr 1 = r 1 K = i, [E,F] = K ^ K ^ , (A.l) 

q — q 

KEK- 1 =q 2 E, KF JT 1 = q~ 2 F . (A.2) 

It is a Hopf algebra deformation of the universal enveloping algebra of s?2 with the comulti- 
plication 

AK = K®K, AE = E^l + r 1 ®E, AF = F ® K + 1 ® F . (A3) 

Iterating this comultiplication, and evaluating the algebra in its two dimensional repre- 
sentation gives the action of U^slz) on the space of states of a spin chain of length N with 
spin s = |: 

F = E f 3 ® " ' ' ® f 3 ® a ~ ® 1 ® ' ' ' ® 1 = E L v a y ' ( A - 4 ) 
y v ^ ' y 

y- 1 

E = ^1® ■ ■ ■ ®t®cr + ® q-" 3 ®---® q' * , (A.5) 



V N- y 

K = q ff3 ® ■ ■ ■ ® q a3 , (A.6) 
v v ' 

N 

where o~ ± , o 3 are the Pauli matrices 

^(oJ)< —(!!!). '-(J- ,). ^ 



1°>=il „-l|- (A.S) 



and 

n a3 = ( 

v q 

In our case, q is not a root of unity, and therefore (C 2 ) 0N splits into the direct sum of 
irreducible highest weight representations of ii^sfc)- 

B The Bethe ansatz for the finite XXZ chain 

In this appendix we recall the expressions for the Bethe vectors for the finite chain of length 
N with UAslj) -invariant boundary conditions. 

Special boundary conditions for the XXZ Hamiltonian when it commutes with the action 



I A. 4 1 of Uq(sl2) on the space of states were found in (7|. Even more remarkable is that this 
Hamiltonian describes an integrable system: the Bethe alsatz for it was constructed in J7|, 
the algebraic version and the complete systems of commuting integrals were found in [9|. 
The Hamiltonian found in [7J has the following form: 

1 N/2-1/2 

w = - o E KV+i + °k°l + i + a + \)i4 - 4+i) - a) • (b.i) 

k=-N /2+1/2 
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Eigenvectors for this Hamiltonian in the subspace where m spins are down and the rest 
up are given by the following formulae: 



B( N )(z) 



£ AW(x|z)Q( N )(x) 



(B.2) 



where Q( N ' (x) is a state with m spins down in the positions with coordinates x, and 

aW(x|z) = £ a (N) (^ ^,J4\ •••4:, (b.3) 



mew 



Here the sum runs over the permutations and reflections (z,- i— > z- x ) of { Zi, . . . , z m }, and 
is the signature of w. The coefficients (z) are: 



where 



and 



;'=1 l<;'<Km 



B(z 1/ z 2 ) = (1 -2Az 2 + Z!Z 2 ) (l -2AZ" 1 +Z2Z" 1 ) , 



j 6(z)= 1-^z-Mz 



1 1 r (N+l)/2 



(B.4) 



(B.5) 



(B.6) 



The vector (B.2 1 is an eigenvector if the numbers Zv have modulus Z/ = 1 and satisfy the 
Bethe equations: 

j?LB(zr 1 /Z/ ) 



,2N 



n 



1=1 B(zj,z t ) 



-, j = l,...,m. 



(B.7) 



The corresponding eigenvalue is given by 



fW(*) = £ 2A-z ; -zri . 

7=1 



(B.8) 



It is quite remarkable that in addition to being eigenvectors these vectors are also Uq^slz) 
highest weight vectors with the weight N — 2m [7|: 



£B( N )(z) =0, 



KB^ N \z) = q N - 2m B^(z). 



(B.9) 



Since the Hamiltonian (B.l i commutes with the action of Uqislz), the vectors FPB^ N \z) are 
also eigenvectors. 
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